gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_10.m

    s=tf('s');
Gk=1/s/(0.5*s+1)/(s+1);
G=feedback(Gk,1);
w=logspace(-1,1);
[mag,phase,w]=bode(G,w);
[Mp,k]=max(mag);
resonantPeak=20*log10(Mp)
resonantFreq=w(k)
n=1;
while 20*log10(mag(n))>=-3
n=n+1;
end
bandwidth=w(n)
bode(G,w),grid;